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We describe a hybrid laser-microwave scheme to implement two-qubit geometric phase gates in crystals of 
trapped ions. The proposed gates can attain errors below the fault-tolerance threshold in the presence of thermal, 
dephasing, laser-phase, and microwave-intensity noise. Moreover, our proposal is technically less demanding 
than previous schemes, since it does not require a laser arrangement with interferometric stability. The laser 
beams are tuned close to a single vibrational sideband to entangle the qubits, while strong microwave drivings 
provide the geometric character to the gate, and thus protect the qubits from these different sources of noise. A 
thorough analytic and numerical study of the performance of these gates in realistic noisy regimes is presented. 
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I. INTRODUCTION 

Quantum information processing (QIP) holds the promise 
of solving certain computational tasks more efficiently than 
any classical device [1]. This prospect has stimulated an enor- 
mous technological effort, whereby prototype quantum pro- 
cessors based on different technologies have already been de- 
veloped [2]. The promised supremacy of these processors re- 
lies on the vast parallelism available at the quantum realm, 
which in turn rests on the quantum superposition principle. 
Therefore, any quantum processor must be well-isolated from 
environmental sources of noise, since these tend to degrade 
quantum superpositions through the phenomenon of decoher- 
ence. Additionally, to benefit from the large parallelism, quan- 
tum processors should also allow for a very accurate manipu- 
lation of the information. Experimental imperfections in such 
manipulation will also conspire to reduce the potential of QIP. 

Hence, the design of QIP protocols that are robust to the 
most important sources of environmental and experimental 
noise, is considered to be a task of primary importance. The 
particular level of robustness required, which can be quan- 
tified by the maximal allowed error e per step of the rou- 
tine, is determined by the possibility of implementing proto- 
cols of quantum error correction that allow for fault-tolerant 
QIP [3]. The so-called fault-tolerance threshold, Spj, typi- 
cally depends on the dominant source of noise, the particu- 
lar error-correcting scheme, technological limitations of the 
experimental platform, and the particular purpose of the rou- 
tine (i.e. initialization, manipulating information, or measure- 
ment). The manipulation of the information is argued to hold 
the most stringent thresholds [4], and although some error- 
correcting schemes with thresholds as high as £fT[ ~ 10~ 2 ex- 
ist, it is commonly agreed that reducing errors below £ft 2 ~ 
10~ 4 is a guiding principle for the development of quantum 
processors. Since the information is usually stored in the so- 
called qubits (i.e. two-level systems), these routines can be 
divided in prescribed sequences of one- and two-qubit gates 
allowing for universal quantum computation [1]. Therefore, 
the quest is to reduce the error of both types of gates below 
the fault-tolerance threshold in the presence of noise. 

Among the different existing platforms for QIP, small crys- 
tals of trapped atomic ions in radio-frequency traps [5] are 
considered to be one of the most promising devices where 



these demanding thresholds could be achieved in the future. In 
fact, one-qubit gates with errors £i q < £ft 2 have already been 
demonstrated using trapped-ion hyperfine qubits subjected to 
microwave radiation [6]. Additionally, two-qubit gates with 
errors as low as £ft, < £2 q < £ft! have also been demonstrated 
for trapped-ion optical qubits subjected to laser radiation [7]. 
Since both qubits operate in a different frequency range (i.e. 
microwave/radio-frequencies versus optical), the reached ac- 
curacies have not been combined in a single device yet. Be- 
sides, it would be desirable to lower the error of the two-qubit 
gates even further to achieve £2 q < £ft 2 • 

In this work, we describe theoretically a scheme capable 
of achieving the two aforementioned goals. Our scheme is 
particularly designed for hyperfine trapped-ion qubits, enjoy- 
ing thus the advantages of using microwave radiation for one- 
qubit gates. In addition, we introduce a scheme for a driven 
geometric phase gate acting on two qubits, which combines 
the advantages of microwave and laser radiation. This hybrid 
laser-microwave scheme can attain errors as low as £2 q < £ft 2 
in the presence of (i) thermal noise, (H) dephasing noise, (Hi) 
phase noise in the lasers, and (iv) intensity fluctuations of 
the microwave driving. In comparison to our previous pro- 
posal [8], the protocol presented in this work allows us to in- 
crease the gate speed by at least one order of magnitude as a 
result of working closer to particular resonances of the laser- 
driven couplings. More importantly, such a near-resonance 
regime turns out to be the parameter regime where the re- 
cent experimental demonstration [9] of the two-qubit gate [8] 
has been performed [10]. Therefore, the present work will 
be useful to provide a theoretical background for the results 
presented in [9], where £2 q ~ £ft! was achieved. Our analyt- 
ical and numerical results show how the geometric character 
of these laser-microwave driven gates underlies the gate re- 
silience to several sources of noise. Moreover, we optimize 
the laser-microwave parameters, such that the errors can be 
reduced below the fault-tolerance threshold £2 q < £ft 2 - 

This article is organized as follows: In Sec. II, we de- 
scribe the mechanism for the driven geometric phase gates 
by means of detailed analytical analysis supported by numer- 
ical results. The relevant analytical calculations are presented 
in Appendix A. In Sec. Ill, the scheme is tested against sev- 
eral sources of noise (i.e. thermal, dephasing, phase, and in- 
tensity noise). Some of these noise sources can be modeled 
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as a stochastic Hamiltonian term, whose properties are de- 
scribed in Appendix B. Finally, we present some conclusions 
and prospects in Sec. IV. 

H. DRIVEN SINGLE-SIDEBAND GEOMETRIC PHASE 
GATES 

A. Two-ion crystals as the hardware for quantum logic gates 

Let us start by describing the system under consideration: 
a two-ion (N = 2) crystal confined in a linear Paul trap [11]. 
Under certain conditions [12], such radio-frequency traps pro- 
vide an effective quadratic confining potential, which is char- 
acterized by the so-called axial co z , and radial {co x ,(Oy} trap 
frequencies. Moreover, when {(0 x ,(0y} (0 Z , the ion equi- 
librium positions arrange in a string along the trap z-axis. 
As customary in these cases, when the particles only per- 
form small excursions around their equilibrium positions, 
their motion can be described in terms of collective vibra- 
tional modes, whose quantized excitations lead to the well- 
known phonons [13]. In the case of ions, the low-temperatures 
provided by laser-cooling techniques justify such a treat- 
ment [14], and the vibrations of the crystal are described by 

N 
a n =\ 

where aj a (a„ a ) are the creation (annihilation) phonon op- 
erators, (Q a , n is the vibrational frequency of the n-th normal 
mode along the axis a = {x,y,z}, and we set h — 1 through- 
out this manuscript. According to Eq. (1), the vibrations along 
different axes are decoupled, which allows us to focus our at- 
tention on the radial modes along the x-axis. However, we 
emphasize that our proposal also holds for the other modes. 
To simplify notation, we write a nx — > a„, co„. x — ► (O n . Note 
that for a two ion crystal there are two radial normal modes: 
the center-of-mass (com) and the zigzag (zz) mode which we 
also shall denote modes "1" and "2", respectively. 

In addition to the vibrational degrees of freedom, our two- 
ion crystal also has internal (i.e. atomic) degrees of freedom. 
We consider ion species with a hyperfine structure, such that 
we can select two levels from the ground-state manifold to 
form our qubit {|0,), This particular choice of qubit has 
two important properties: (i) spontaneous emission between 
the qubit states is negligible (i.e. T\ times are much larger than 
experimental time-scales), and (ii) the typical qubit frequen- 
cies <x>o/2k = (E\ —Eq) /2k ^1-10GHz allow for microwave, 
or radio-frequency, radiation to drive directly the qubit tran- 
sition. To test our scheme with experimental realistic param- 
eters, we shall consider 25 Mg + in this work, although we re- 
mark that our scheme also works for other ion species (see 
Fig. 1). We choose two hyperfine states |0) = \F = 3,M = 3) 
and |0) = \F' = 2,M' = 2), where M,M' are the Zeeman sub- 
levels of the F, F' levels of the electronic ground state mani- 
fold, which are separated by a frequency (Oo/2n « 1.8 GHz. 

By considering the standard Zeeman interaction between 
the atomic magnetic-dipole and an oscillating magnetic field 



from the microwave source, and using the standard rotating- 
wave approximation in quantum optics [15], the qubit Hamil- 
tonian can be written as follows 

N 

H q = £|fflbof + ^(n d o^e- ifl *i'+Hx.), (2) 

where we have introduced the qubit operators af = — 
|0;)(0,-|, ct+ = |1 ; )(0,| = {a^)\ and the frequency (Rabi fre- 
quency) 0)d(£id) of the microwave driving. Note that the ro- 
tating wave approximation requires w Oo, and \Q.d\ <C COq, 
which still leaves the possibility of achieving very strong 
drivings Q.^/271 ^1-10 MHz, while benefiting from the inten- 
sity stability of microwave sources. We also note that the 
large wavelengths of the microwave traveling waves avoid 
any qubit-phonon coupling induced by the microwave, which 
would impose further constraints on the strengths of the driv- 
ings. In this work, we exploit the possibility of obtaining such 
strong microwave drivings to produce fast two-qubit geomet- 
ric phase gates with robustness to several sources of noise. 

As done in previous two-qubit gate schemes [16-22], the 
main idea is to consider an additional coupling between the 
qubits and the phonons, and use the latter to mediate inter- 
actions between distant qubits in the ion trap. For the spe- 
cific qubit choice, such qubit-phonon couplings can also be 
achieved in the near-field of microwave sources, by produc- 
ing either static [23, 24], or oscillating magnetic-field gra- 
dients [25, 26]. The generality of the scheme presented 
in the following subsection would allow us to use any of 
these approaches for an all-microwave protocol. However, 
we follow [8], and consider laser-induced qubit-phonon cou- 
plings [11], which have already been demonstrated in several 
laboratories. Therefore, our work can be considered as an in- 
stance of a hybrid laser-microwave protocol. 

The qubit-phonon coupling is provided by a two-photon 
stimulated Raman transition via a third auxiliary level. The 
transitions to the auxiliary excited state are driven by two far 
detuned laser beams in a traveling-wave configuration, such 
that their detuning is much larger than the excited-state de- 
cay rate (e.g. for 25 Mg + the linewidth of the excited state 
is r/27T « 4 1.4 MHz, and thus a detuning on the order of 
A/2n « 10-100GHz would suffice, see Fig. 1). In this limit, 
the excited state can be adiabatically eliminated [27], and the 
qubit-phonon coupling becomes 

tfqp = L ^a+c^-^) +H.c, (3) 

i 

where Hl = Q.l 1 Q.l 2 /2A is the two-photon laser Rabi fre- 
quency, and the two-photon wavevector (frequency) kL = 
kL, — kL 2 ((Ol = ©Lj — C0l 2 ) stre defined in terms of the corre- 
sponding parameters of the two laser beams [28]. By directing 
kL along the x-axis, and setting the laser frequencies such that 
0>l ~ COq — (On, and the detunings 8 n = (0^ — (coq — CO„), such 
that 8„ <C (0„ and | £ir_ | -C (0„, a rotating-wave approximation 
leads to the so-called first red-sideband excitation 

ffqp = ^^ in ^a n e-^' +H.c. (4) 

i.n 
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Figure 1 : Ingredients for the driven geometric phase gate: For a generic ion species with hyperfine structure, as determined by the quantum 
numbers n 2 LjF, the qubit is chosen from the two hyperfine states F = /± j of the ground-state manifold n 2 S\n- Additionally, a quantizing 

magnetic field splits the F-states into Zeeman sub-levels, such that |0) = \F =1+ \ ,M) and |1) = \F' = 1 — \ ,M') define the hyperfine 
qubit. Two laser beams with Rabi frequencies Q.^ , Q.^ drive the first red-sideband transition of the qubit in a stimulated Raman configuration 
|0)c3>|tt) — > |1) §3 \n— 1), where n is the number of phonons of a particular mode of the two-ion crystal. Simultaneously, the qubit transition is 
driven directly by a microwave 



Here, we have defined the red-sideband coupling strengths 



of producing two Bell states at particular instants of time 



2V2 



^■2 = i(-l) 



|fl L |T72 

2V2 ' 



(5) 



where we have set the Raman-beam phase <f>L = 0, and 
used the Lamb-Dicke parameters 77,, = kL • e. v / \/2m(Q„ <C 1. 
We will refer to these couplings generically as forces, since 
^i n \J2m(S) n corresponds to a force applied on the harmonic 
oscillator representing the vibrational mode. 

As explained below, equations (1), (2) and (4) form the 
driven single-sideband Hamiltonian 



Hds; 



Hp+H q +H qp7 



(6) 



which can be used to obtain the desired geometric phase gates 
for strong-enough microwave drivings. At this point, it is 
worth commenting that strong qubit drivings have also been 
realized experimentally in combination with a state-dependent 
force [29]. The driving in this case endows the scheme with 
protection from dephasing noise, while the resilience to the 
thermal motion of the ions is provided by the state-dependent 
force. However, this scheme is prone to phase noise due to 
fluctuations in the laser-beam paths, or to intensity fluctua- 
tions of the qubit driving. In a different context, we should 
also mention that strong-driving-assisted protocols for the 
generation of entanglement have also been considered theoret- 
ically for atoms in cavities [30], where the driving may reduce 
errors due to thermal population of the cavity modes [31]. 



B. Driven geometric phase gates 

1. Introduction: geometric phase gates with a single sideband 

In reference [8], we have analyzed how the driven single- 
sideband Hamiltonian H& %% leads to an entangling gate capable 



|01)^|¥-> = i(|01)-i|10>), 
|10)^|¥+) = i(|01)+i|10}), 



(7) 



in the far-detuned regime \&i n \ <§C 8„. In this limit, the 
weak qubit-phonon coupling (4) is responsible for second- 
order processes where phonons are virtually created and an- 
nihilated. This term alone leads to a flip-flop qubit-qubit 
interaction i/j nt °= Jf^a^aJ +Hc. via virtual phonon ex- 
change, where the coupling strength is given by Jfj 1 = 
— J?,„i^*,/(S„ with the above introduced forces and detun- 
ings jF !n and 8„, respectively. At certain instants of time this 
interaction can act as an entangling gate. However, we also 
showed there that the performance of such a gate would be 
severely limited by thermal and dephasing noise. The main 
idea put forward in [8] was to exploit a strong resonant mi- 
crowave driving (2), fulfilling cod — coq and 8„ <C ild £ K, 
as a continuous version [32, 33] of refocusing spin-echo se- 
quences [34, 35], to protect the two-qubit gate from these 
sources of noise. In this case, the virtual phonon exchange 
leads to an Ising-type interaction H mt °< Jf^afaj + H.c, 
where we have introduced of = + af and If? — 1 /47??. 
This qubit-qubit interaction can generate all of the four Bell 
states 



100} -»• I*") = ^(|00> -isgn(/ff)|ll>), 
|01>^| l P-> = ^(|01>-isgn(/^)|10)), 
|10) -> |^+> = ^(|01) +isgn(./] 2 ff )|10)), 
|H) ^l* + > = 75(100) +isgn(y]f)|ll)). 



(8) 



Unfortunately, the explored far-detuned regime leads to gates 
that are more than one order of magnitude slower than 
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state-of-the-art implementations based on other schemes [36]. 
Therefore, although both the simplicity of the gate, and its re- 
silience to different sources noise are interesting advantages, 
its lower speed presents a considerable drawback. 

In this work, we show how one can abandon the far-detuned 
regime, while still preserving the nice properties of the driven 
single-sideband gate. Below, we show that by working in the 
context of the geometric phase gates [18, 20, 37], our driven 
nearly -re sonant single-sideband gate: ( i) Can attain speeds of 
tg ~ 63 jits, which is one order of magnitude faster than the far- 
detuned gate [8]. (ii) Can minimize thermal and dephasing er- 
rors down to £ t h, £d < 10~ 4 , which improves by more than one 
order of magnitude the far-detuned gate [8]. (Hi) Can with- 
stand fluctuations of the laser phase occurring on timescales 
longer than the 63 jxs, such that e p h < 1CT 4 , which was out- 
lined in [8], but not analyzed carefully, (iv) Can also resist 
relative fluctuation in the intensity of the microwave driving 
at the 10~ 4 -level directly with errors E\ ~ 10~ 3 . Moreover, we 
show that by adding a secondary driving, the scheme can sup- 
port stronger intensity fluctuations, while providing smaller 
gate errors £i < 10~ 4 . 



2. Qualitative analysis: dressed-state interaction picture and 
rotating-wave approximation 

To understand the mechanism underlying the driven nearly- 
resonant single-sideband gate, let us make a small detour, 
and consider the so-called geometric phase gates by state- 
dependent forces [18, 20, 37]. By combining the red-sideband 
(ofa„) term (4) with a blue-sideband (<7+aJ) that has an 
opposite detuning but equal strength and adjusting the laser 
phases appropriately [18], the qubit-phonon Hamiltonian be- 



tangled, such that the evolution operator becomes 



" 8 n 



J" 



(11) 



This unitary can be easily seen to provide the desired entan 
gled states (8) when ^(27^) = n/4. Additionally, if the initia 
spin states are eigenstates of G^G^ , the gate gives the table 



-e'2 



(12) 



up to an irrelevant global phase. This corresponds to a two- 
qubit 7r/2-phase gate which, together with single-qubit rota- 
tions, gives a universal set of gates for quantum computation. 
The fact that the phase-space trajectory is closed at tg, allows 
for the interpretation of the 7r/2-phases as geometric Berry 
phases determined by the area enclosed by the trajectory [37]. 
Let us remark two properties of these gates: (i) They do not 
rely on any far-detuned condition | J^,„| <C 8„, and can be thus 
much faster, (ii) The geometric origin of the phase gate under- 
lies its robustness with respect to thermal motion of the ions 
(i.e. spin-phonon disentanglement occurs when the trajectory 
closes regardless of the vibrational state). 

After this small detour, it becomes easier to understand why 
a gate based only on a nearly-resonant single-sideband cannot 
lead to a geometric phase gate, and is thus very sensitive to 
thermal noise. By rewriting the single-sideband (4) in the in- 
teraction picture with respect to Hq for Q.^ = 0, we get 



qp 



I ar. 

-> it 



[of + iof)a n e- iS "' +H.c. 



(13) 



qp ' 



a„e 



-i8„t 



He, 



(9) 



where the "hat" refers to the interaction picture with respect 
Ho = H q + Hp, and we have switched off the microwave driv- 
ing = 0. This hamiltonian (9) can be understood as a push- 
ing force in a direction that depends on the qubit state in the 
x-basis of\±i) x — ±|±,) x . Since a single Pauli matrix appears 
in the above equation, the time evolution under such a state- 
dependent pushing force reduces to that of a forced quantum 
harmonic oscillator [38], which can be solved exactly (see 
e.g. [39]). This leads to the following time evolution opera- 
tor 



U{t) 



a E,„(^f(e-*'-l)afa„-H.c.) 



jn i f sin(5„>) 



(10) 

The first unitary corresponds to a displacement operator, 
which is well-known in quantum optics [15], with the pecu- 
liarity of being state-dependent. In phase space, this term in- 
duces periodic circular trajectories for the vibrational modes 
that depend on the collective spin state, and thus leads to 
qubit-phonon entanglement. When f g = k n 2n/8 n , k n E Z, the 
trajectories close and the qubits and phonons become disen- 



where a] = — i(7 ( - +i<7^. Hence, it is clear that the red- 
sideband yields a combination of two state-dependent forces 
acting on orthogonal bases. This fact forbids an exact solu- 
tion, such as the one obtained for a single state-dependent 
force (10), and we inevitably lose the notion of state- 
dependent trajectories that close independently of the vibra- 
tional state. Accordingly, the qubit and phonons get more en- 
tangled as the far-detuned condition <C 8 n is abandoned, 
and the gate fidelity drops severely for ions in thermal motion. 

If we now switch on a resonant microwave driving ©d = 
O)o, fid € K, then the Hamiltonian (13) changes into 

H qp=L^r( (7 -+ i ^ cos(n d f)-i(7fsin(n d f))fl„e- iS "'+H.c.. 

i.n l 

(14) 

Here, the "tilde" refers to the following "dressed-state" inter- 
action picture, namely //qp = \j&(t)H 0S lj'^{t), where 



It is now possible to argue that the state-dependent a y and a z 
forces rotate very fast for sufficiently strong drivings ^> 8 n , 
and can be thus neglected in a rotating wave approximation 
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provided that |^,„| <C £2 d - Under this constraint, the driven 
single-sideband Hamiltonian becomes 



qp ' 



^{■^inOfant-^' 



H.c. 



(16) 



and thus a formal solution like Eq. (10) becomes possible 
again, provided that we make the substitution & 



in 7 o ln ' 

According to this qualitative argument, a strongly-driven 
single-sideband Hamiltonian can also produce the desired ge- 
ometric phase gate. Therefore, it seems possible to abandon 
the far-detuned regime |J^,„| <C 8 n of our previous work [8], 
and obtain faster gates that are still robust with respect to ther- 
mal noise. In the following sections, we will present a quanti- 
tative detailed analysis to test the validity of this idea. 



3. Quantitative analysis: Magnus expansion and numerical 
analysis 

As discussed above, the presence of different state- 
dependent forces in the Hamiltonian (13) avoids an exact so- 
lution for the unitary time evolution operator. However, since 
we are interested in the strong-driving limit |J£,„| <C £i d , we 
can use the so-called Magnus expansion [40], truncating it to 
the desired order to obtain the leading contributions to the dy- 
namics (see Appendix A for the details). The time-evolution 
operator in the Schrodinger picture reads 



f/app(0 = f/ d t (0e a(0 , « a, (/)+ o 2 (/ ) + z), 



and ^ = (i2Li? n )"/^d' X — (^L^J«) /^d$n are tne small pa 
rameters in the Magnus expansion. In this expression, we 
have defined the anti-unitary operators containing the differ- 
ent state-dependent displacements 

^iW=I§ 1 (e- iV -l)afa„-H.c. 



■ S 11 



f£ 4(£2 d - 8„ 

- 7 III 



_ (e .(n d -5„)r _ i)(-iaf + af)a„-H.c. 



(18) 



and the operators leading to the qubit-qubit interactions 

O 2 (0 =i£ 



^■^jn ,,f. sin 8„t 
45„ ff ' ff U r 8„ 



-it £A£2,„ \a\a n - i j of + £ (/„„,(?) of a^a„ - H.c), 

(19) 



where we have introduced the couplings strengths 



Ail,, 



Win? ( 1 



1 



4 \Q. d -8„ i2 d + <5„/' 



(20) 



and the following time-dependent functions 



fnm(t) : 



8(5„-5 ffl )Vn d - 



-V + ^-^Ve^")'-!). 

-8 m n d + 5 m y v > 



(21) 

At this point, it is worth comparing the final expres- 
sions (18)-(19) to those corresponding to a state-dependent 
force in Eq. (10). The first-order contribution, (18), resem- 
bles the state-dependent displacement in (10), but we get in 
addition the contribution of more state-dependent forces in 
different bases. From the first line of the second-order contri- 
bution (19), we observe also a similarity with the qubit-qubit 
couplings in (10), but we get additional residual qubit-phonon 
couplings in the second line of (19). We discuss below how 
all these additional terms can be minimized, such that we are 
left with an effective time evolution that is analogous to that 
of a state-dependent force. 

Let us note that Eqs. ( 1 8) and (19) correspond to the leading 
terms in the second-order expansion expansion of U(t ), while 
the complete expression may be found in Appendix A. Before 
analyzing how the geometric phase gates arise from Eqs. (18)- 
(19), let us check numerically the validity of our derivation by 
comparing it to the time-evolution under the full driven single- 
sideband Hamiltonian // dss (6). To perform the numerics more 
efficiently, we expressed // dss in a picture where it becomes 
time independent, namely 



£ §na\a n ■ 



E 



CTf + ^(^ m a+a n + He), (22) 



(I 7 ) where H' dss = U' (t)H iss (U' (t))\ and the unitary is 



U'(t) — e i '£"( <H o~ (B L) a " a "e i ' £ ' i* 00 ®! 



(23) 



Accordingly, the time-evolution operator can be written as 



f/exact(f) = (^(0) t e^'. 



(24) 



For the numerical simulations, we chose realistic parameters 
for ion-trap experiments, which are summarized in Table 1. 
The outcome of the simulations is shown in Fig. 2(a), which 
shows a very good agreement between the Magnus expansion 
and the exact time evolution for the qubit dynamics. This sup- 
ports the validity of our derivations, and allows us to carry on 
with the description of the driven geometric phase gate. 

The action of the first contribution i2i(f) (18) to the time- 
evolution operator (17) can be understood as follows: for 
each vibrational mode, three non-commuting state-dependent 
forces aim at displacing the ions along circular paths in phase 
space, such that the specific direction of each trajectory de- 
pends on the particular eigenstates |±. v ), |± v ) and |0/1) of the 
operators cr x , a y and O z . For weak drivings, the combina- 
tion of these forces will deform the circular paths, and lead 
to phase space trajectories that are not closed any longer (see 
Fig. 3(a)). It is clear from the analytical expression (18) that, 
by applying a sufficiently-strong driving i2 d ^> 8„, \^j n \, the 
forces in the a y ' z -basis get suppressed, and, to a good approx- 
imation, we are left with the desired single state-dependent 
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Table I: Values of trapped-ion setup for the numerical simulation 

cq : /2x a> x /2n 8 mm /2it r) com 8 zz /27t r\ zz Q. L /2it Q. A /2n 
1MHz 4MHz 127kHz 0.225 254kHz 0.229 811kHz 7.2 MHz 



three state-dependent displacements in (18), and the residual 
qubit-phonon entanglement introduced in the second term of 
the second line of Eq. (19). In order to achieve this, we need to 
adjust the Hamiltonian parameters such that, at the gate time 
f g , the following constraints are fulfilled 

2% 

t g = r^, 82=k8 u Q. d =p8i, r,k,pe1. (25) 

In other words, we need the detuning of the zigzag mode 82 
and the microwave Rabi frequency Q. d to be an integer multi- 
ple of the center-of-mass mode detuning 8\ . For these condi- 
tions, i2i (f g ) = vanishes exactly. Moreover, we reduce con- 
siderably the qubit-phonon entanglement due to the second- 
order term £22(fg) = — ifg^dss, where we have introduced the 
Hamiltonian 

i,j i,n 
and the following coupling strengths 

n HU « 

Note that a closer inspection of all the additional terms in the 
Magnus expansion (see Appendix A) shows that they all can- 
cel exactly under these constraints. Therefore, all the dynam- 
ics for the gate is exactly described by the Hamiltonian (26). 

Finally, in the limit of very strong drivings £2d ^ 8„, the 
residual qubit-phonon coupling is minimized as AClj n (20) de- 
creases with increasing driving. Moreover, one can even can- 
cel exactly the contribution of the second term in Eq. (26) by 
introducing a spin-echo pulse [8] or, equivalently, a phase re- 
versal [9] of the microwave driving at half the expected gate 
time. As an additional benefit, this procedure guarantees the 
closure of phase space trajectories in the case of an imperfect 
detuning if the forces in the G x ' y bases can be neglected in a 
rotating wave approximation. This has been demonstrated in 
the recent experiment [9]. Therefore, the full time evolution 
can be written approximately as 

U m (h) = Ul(t g )e- [H ^, H dss = Y^jfrcricrl (28) 

y 

From this expression, we can extract the gate time to be 
(27j 1 2 S ) f g = which will give rise to the entangling ta- 
ble (8), or to the phase-gate table (12), depending on the cho- 
sen basis for the initial states. In Fig. 2(b), we check the valid- 
ity of this description by computing the fidelity of generating 
the Bell state \ X P + ) at the expected gate time, starting from the 
initial qubit state 1 10). As seen from the figure, the agreement 
between the full Magnus expansion and the exact dynamics is 
very good, and we obtain the desired entangled state at f g with 
very high fidelities. 

Note that we have not made any assumption on the ratio be- 
tween the sideband strengths and the laser detunings | J^„ |/6"„, 
and thus we do not require to work in the far-detuned regime 
\^in\/8 n <§C 1 of our previous proposal [8]. According to the 
expression of the gate coupling strengths jfj s , this means that 




t{lls) 120 




tQls) 120 

Figure 2: Entangling ofcr? gate: The figure displays the exact 
dynamics by numerical integration of the unitary evolution (24) 
(solid lines) with the dynamics given by the approximated time- 
evolution operator (crosses). We consider an initial state po = 
|10)(10| ®Pth("li"2), where hi ~ «2 ~ are the initial mean num- 
ber of phonons for the cm and zz modes. Note that due to the nearly- 
resonant sideband, phonons are created during the gate, and we have 
to set a high truncation n max = 10 to the vibrational Hilbert spaces. In 
(a), we represent the dynamics of local qubit operators (cr?), which 
show the qubit flip 1 10) -> |01) after t = 126 /is. In (b), we dis- 
play the fidelity between the time-evolved state and the Bell state 
|<P + ) = (|01) +i|10))/\/2. At t = 63 /is, this fidelity approaches 
unity. The agreement of both descriptions supports the validity of 
our analytical derivation. 



force in the 0" T -basis with a halved strength J£",„ — > \-^in with 
respect to the standard M0lmer-S0rensen term (9). Accord- 
ingly, we should recover the circular phase-space trajectories 
when the initial state of the qubits is an eigenstate of the a\ oi 
operator. In Fig. 3(b), we analyze the strong-driving dynam- 
ics numerically, and confirm the above prediction. Therefore, 
the qualitative description of the previous section is put on a 
firmer ground by the use of the Magnus expansion. 

The Magnus expansion also allows us to improve further 
the geometric character of the gate. By setting the correct 
values to the detunings and intensities, it is possible to can- 
cel exactly the qubit-phonon entanglement introduced by all 
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that 82 = kS\ with k E Z, we get 



g-1 



(30) 



where | = ©2/01 < 1 is the ratio of the normal-mode frequen- 
cies. To make the gate as fast as possible f g = 2tz/(\c\coi), 
the modulus of c should be as big as possible, and we thus 
choose k = 2, so that |c| = 1 — ^. According to these con- 
siderations, the gate time in units of the trap frequency is 
fgOjc = 2n/(\ — ^ ), which is minimized by minimizing B, (i.e. 
large difference of the center-of-mass and zigzag mode). This 
occurs when approaching the linear-to-zigzag instability. 

Finally, the laser intensities in jff nave to be adjusted such 
that we can fulfill the following condition 



2n 



(31) 



-0.4 (i cm > 0.4 



where r £ Z is an integer that determines how many loops in 
phase space the center-of-mass mode performs. On the one 
hand, r should be as small as possible to maximize the gate 
speed. On the other hand, it turns out that the smaller r is, the 
more fragile the gate becomes with respect to dephasing noise 
(see sections below). Therefore, we optimize the value of r to 
find a compromise between gate speed and gate fidelity. 
Inserting expression (27) in the above equation, we obtain 



Figure 3: Phase-space trajectories: Numerical analysis of the 
phase-space trajectories under the full driven single-sideband Hamil- 
tonian (6). (a) Trajectories of the forced normal modes in the weak 
driving regime 0.^/2:1 = 200kHz. In the left panel, the initial spin 
state is | + +) or | ) while the phonons initially are in a ther- 
mal state Pth^l,'^) = Pth(«l) ® Pth(«2) with hi « h\ = 0.5, and 
according to the forces (5), only the center-of-mass mode is driven 
in phase space (x cm ) = 1/ \p2(a\ +«}), (Pcm) = i/\ // 2{a[ — a\). In 
the right panel, the initial spin state is | H — ) or | — h), such that only 
the zigzag mode develops a driven trajectory in this case (x zz ) = 

l/\/2(a2 +^2)1 (.Pzz) = V 1 ^-{ a \ ~ a2 )- In both cases, the trajecto- 
ries are not closed, (b) Same as (a) but in the strong-driving limit 
£\)/27T = 5 MHz. As an effect of the strong driving, the trajectories 
are closed, and the gate inherits a geometric character. In the numer- 
ics we have to set a truncation rc max = 15 to the vibrational Hilbert 
spaces. 



we can increase the gate speed considerably by working closer 
to the resonance (i.e. f 2 w 63 /is). The conditions (25) also im- 
ply that the phase space trajectories of the vibrational modes 
are closed, providing a geometric character to the two-qubit 
gate, and its resilience to the thermal noise. 

In the following sections, we will analyze the speed 
and noise robustness of this driven nearly-resonant single- 
sideband gate. However, let us first address how the con- 
straints (25) can be met in practice. The detunings for the 
center-of-mass and the zigzag modes are given by 

81 = col- (fflo-ffli), 8 2 = (Ol- (0J0-OJ2). (29) 

By setting 8\ = cco\ , such that c is a constant, and demanding 



Hi 



(32) 



which fixes the laser intensities. Finally, we note that the mi- 
crowave driving strength can be modified in the experiment 



&d = p8u 



(33) 



with p 6 Z being a very large integer to meet the strong- 
driving condition Q.^ 3> 8„. With these expressions, the pa- 
rameters of the setup are fixed in terms of two integers (r,p). 
In the following section, we optimize the choice of these inte- 
gers to maximize simultaneously the gate fidelity and speed. 
We explore how such fidelities approach the fault-tolerance 
threshold in the presence of different sources of noise as the 
microwave driving strength is increased. 



III. GATE ROBUSTNESS AGAINST DIFFERENT 
SOURCES OF NOISE 

In this section, we provide a detailed study of the behavior 
of the gate in the presence of noise. We will analyze four dif- 
ferent sources of noise: (i) Thermal noise, which is caused 
by the thermal motion of the ions (i.e. we do not assume 
perfect ground-state cooling of the radial modes), (ii) De- 
phasing noise, which is caused by fluctuating Zeeman shifts 
due to non-shielded magnetic fields, or by fluctuating ac-Stark 
shifts due to non-compensated energy shifts caused by fluctu- 
ating laser intensities. (Hi) Phase noise, which is caused by 
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fluctuations in the phases of the laser beams, (iv) Intensity 
noise, which is caused by fluctuations in the microwave driv- 
ing strength. We will show that the driven single-sideband 
gate is intrinsically robust to noise of the type (i), (ii), (Hi), 
whereas it will suffer from (iv) at very strong drivings. We 
propose to incorporate an additional weaker driving to make 
the gate also robust with respect to intensity fluctuations. To 
show all these properties, we integrate numerically the dy- 
namics (24) given by the full Hamiltonian (22) subjected to 
the additional noise sources. In the simulations, the parameter 
r in (31) and (32) is set to r = 8, since we have found that it 
gives the best compromise between fidelity and speed. 



A. Resilience to thermal noise 




10 _ ,„ fxjra .30 40 



One of the primary sources of gate infidelities in early 
trapped-ion proposals [16] was the noise introduced by resid- 
ual thermal motion of the ions after a stage of resolved- 
sideband cooling. In order to consider such thermal noise, 
we consider an initial state po = |10)(10| <S> Pth(«i,«2) for a 
thermal phonon state Pth(«i,«2) = Pth("i) <8 Pth(%) charac- 
terized by the mean phonon numbers {ni,«2} for each mode. 
Note that due to the small difference in frequency of the nor- 
mal modes in x-direction we have h\ ~h~2- According to the 
table (8), the target state after the gate is p(? g ) = | <£> 

p th («i,« 2 ), where = (|01>+i|10»/>/2. 

Based on the Magnus-expansion evolution operator (17), it 
is possible to derive a closed expression for the gate fidelity 

J^ +> =Tr{|^ + )(^ + |®l ph f/ M (fg)pof/M(fg)}. (34) 

in the strong driving limit £2^ 3> 8„ (i.e. neglecting all terms 
that are suppressed with 1 /£2d)- After some algebra, the fi- 
delity can be expressed as a sum of expectation values of the 
displacement operator 

(D(a,,)) = Tr ph {e a » a »- a >"p th («„)} = e -l a "l 2 (""+ s). (35) 

Provided that the gate time fulfills t^Q.^ = 4np where p E Z, 
the fidelity of the entangling operation reads as follows 

=\ + le~Z**'( 1 - aa W> + i£ n e-4*H(i-co S (iyg)) 

(36) 

i or a 

where K n = 1 'g 1 (2n„ + 1 ) . From this expression, it becomes 
clear that the fidelity would decrease exponentially with the 
mean number of thermal phonons, unless the phase-space tra- 
jectories are closed 8 n t g = 2n. Under such condition, the gate 
fidelity is maximized J^+vffg — 2n/8 n ) = 1 in the limit of 
ild — > oo. This expression not only unveils the importance of 
the geometric character of the gate, but will also be useful in 
understanding the effect of other sources of noise. 

In order to check how the fidelity approaches unity as the 
microwave driving £2<j is increased, we study numerically the 
exact time evolution (24). To avoid timing errors due to the 
fast oscillations induced by the strong microwave driving, we 
add the above mentioned refocusing spin-echo pulse at half 
the expected gate time, such that the total time evolution is 



Figure 4: Thermal noise: Error e = 1 — & of producing the Bell 
state | I F + ) from the initial state 1 10) at the expected gate time t„ as 
a function of the applied microwave driving field fljj. Note that each 
point of the figures fulfills the condition flj = p8\ for an integer 
p. The thermal phonon population of the initial state h\ is increased 
from to 1 . Gate errors well-below the fault-tolerance threshold (FT) 
can be achieved £2 q < £ft 2 . 

Um{tg,G) = f/ e xact(fg, |)(fTiO-|)C/exact(|,0). We then com- 
pute numerically the fidelity of producing the desired Bell 
state 

JV) =Tr{| l P+)('F+|®l ph f/ full (f g ,0)pof/ f l lll (fg,0)}. (37) 

In order to integrate the full Hamiltonian numerically, we 
truncate each vibrational Hilbert space to a maximum num- 
ber of phonons of « max = 25 per mode, which is sufficiently 
high so that no appreciable error is introduced. 

In Fig. 4, we represent the numerical results for the gate er- 
ror e = 1 — & for various microwave driving strengths. By 
setting r = 8 in Eq. (31), we get a gate time of t g s» 63/J.s, 
which improves the gate speed of our previous work [8] by 
one order of magnitude. Besides, the results shown in figure 
4 show that the gate errors surpass the fault-tolerance thresh- 
old £pr, ~ 10~ 4 for sufficiently strong drivings. This figure 
illustrates clearly the benefit of this nearly-resonant single- 
sideband gate in comparison to our previous scheme [8], since 
the gate error even for the largest thermal phonon numbers 
n\ = 1 is reduced by several orders of magnitude. 



B. Resilience to dephasing noise 

Another source of gate infidelities for several trapped-ion 
qubits is the noise caused by fluctuating magnetic fields, or 
laser intensities, which lead to a shift in the resonance fre- 
quency of the qubit through either the Zeeman shift, or an 
uncompensated ac-Stark shift (see Appendix B). On the time 
scales considered, these fluctuations can be modeled as 

tffluc(0=E2- Ac *( f K' (38) 



9 



where A<So(f ) is a normal stochastic Markov process known as 
an Ornstein-Uhlenbeck (O-U) process, and we have assumed 
that the fluctuating source acts globally on the two ions (e.g. 
global fluctuating magnetic field). 

The O-U process, which is defined by a relaxation time z 
and a diffusion constant c [41], is characterized by the follow- 
ing identities 

Acoo(f) =0, Acoo{t)Aa>o(t - s) = f e - *, (39) 

where the "bar" refers to the stochastic average, and we 
assume a Markovian limit where the timescales of interest 
greatly exceed the correlation time. The time-evolution of the 
density matrix in the so-called Born-Markov approximation 

P = - 1°° ds[H &uc (t), [H mc {t - s),p (*)]] , (40) 
Jo 

can be expressed in the following Lindblad form 

p=J? d (p)=r d (S z pS z - l 1 S 2 p-±pS 2 ) 1 T d =^ (41) 

where we have introduced the collective operator S a = L; of 
for a = z. By means of the adjoint master equation [42], one 

finds that of(t) = e^'fff (0) = of (0)e _,CT2 / 2 , which allows 
us to define the decoherence time as T 2 — 2/cT 2 . This ex- 
pression, together with the condition of short noise correlation 
times, fulfilled by setting 1 = 0.1 72, allows us to set the noise- 
model parameters for the study of a variety of decoherence 
times T% £ [15,40]/is. We note that these values are overly 
pessimist, and would only occur for a poor shielding of the 
magnetic fields, or a bad stabilization of the lasers. 

As we show below, even for these poor coherences, the 
strong microwave driving suppresses the noise, and reestab- 
lishes a good coherent behavior. The qualitative argument 
is again to move to the dressed-state interaction picture (15), 
where # fluc (?) = U d (t)H nuc (t)U^(t) becomes 

#fluc(0 = L 5 Afi) o( f ) (cos(iV)af + sin(Q d t)of) . (42) 

i 

For a sufficiently strong driving, the noisy terms rotate very 
fast, and we can neglect them in a rotating-wave approxima- 
tion. 

In order to put this argument on a firmer footing, let us 
use again the Born-Markov approximation (40). We obtain 
a Lindlad-type master equation p = J?? d (p), where the new 
dephasing super-operator is 

^d(p) = ~i[^S x ,p] + ^ £ (S aP S a \S 2 a p ip5 2 J, 

gi=z,v 

(43) 

and we have assumed that Q. d T2 ^ 1. Here, we have defined a 
dressed-state energy shift and a renormalized dephasing rate 

AO — Qd? p _ r d 

Ai2d - 4T 2 (1 + (£2 d T)2) ' d ~ (1 + (£l d T)2) ■ (44) 

The adjoint master equation associated to such Liouvil- 
lian (43) yields the coherences decay of(t) = of (0)e~ 2rit . 
Accordingly, we get a renormalized decoherence time 

f 2 = r 2 (i + (n d T) 2 ), (45) 




Figure 5: Dephasing noise: Error in the generation of the Bell 
state l^ - ) from the initial state 1 00) <g> |«i ,^2) with "« = 0, for a time 
evolution including the dephasing noise. The results shown include a 
statistical average over N a = 10 3 realizations of the noise process per 
point. We set a truncation of the vibrational Hilbert spaces to hold a 
maximum of n max = 7 phonons per mode. The error is shown as a 
function of the applied microwave driving strength, and for different 
values of T 2 £ [15/is,40/is]. Gate errors £2 q below the fault-tolerance 
threshold (FT) can be achieved e 2q < £ft 2 - 



which increases quadratically with the driving strength [43]. 

We now address qualitatively the noise effects on the fi- 
delity of the entangling operation in the strong-driving limit. 
Note that both |10) and lie in the "zero-magnetization" 
subspace, which can be easily checked to be a decoherence- 
free subspace of the original dephasing (41). Therefore, the 
effects of the dephasing noise must be tested for another Bell 
state, such as |<J>") = (|00)-i 1 1 1) ) /a/2, which is gener- 
ated from 1 00) (see Eq. (8)). We identify two possible ef- 
fects: (i) On the one hand, the qubit coherences are degraded 
in a timescale given by T 2 = r 2 (l + (n d T) 2 ). (ii) On the 
other hand, the state-dependent forces for the leading term 
in (17) are damped according to & m — > ^",„e _ '/ r2 (i.e. recall 
of if) = of (0)e~ 2rd '). This avoids the perfect closure of the 
phase-space trajectories, and thus decreases the gate fidelity. 

To validate these predictions, we integrate the full Hamilto- 
nian (22) incorporating the noise term (38) and once again 
a refocusing spin-echo pulse at half the gate time. In Ap- 
pendix B, we describe in detail the conditions under which 
the stochastic noise process can be numerically propagated in 
time, and incorporated into the full time evolution. In Fig. 5, 
we show the numerical results for the gate error as a func- 
tion of the microwave driving strength, where we have set 
h\ = «2 = to distinguish clearly between the effects of ther- 
mal and dephasing noise. This figure shows neatly how the 
dephasing error is considerably suppressed by increasing the 
driving strength. Errors well-below the fault-tolerance thresh- 
old can be achieved again for sufficiently-strong drivings. Ad- 
ditionally, this figure also shows the advantage of the near- 
resonant gates with respect to the far-detuned ones [8], since 
the error is reduced by several orders of magnitude. 
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C. Resilience to phase noise 

Another possible source of noise for trapped-ion gates are 
fluctuations in the laser phases at the position of the ions. 
Such fluctuations become especially dangerous for the geo- 
metric phase gates based on two non co-propagating Raman 
laser beams, which are more sensitive to laser path fluctua- 
tions. This noise can be modeled by substituting the sideband 
couplings ,^i„ — > ,^i n t l ^'\ where A<p(f ) is again a stochastic 
variable. Therefore, the full Hamiltonian (22) becomes 

H' dss = J><4«» +E ^°* + ^ in ^O+a n + H.c). 

n i i,n 

(46) 

However, in this case we model the noise as another nor- 
mal stochastic process, namely a Wiener process A<p <E [0, 2ft]. 
This stochastic process is fully characterized by a single dif- 
fusion constant c and its initial value 

A<p(r) = A<p(0) =0, A<p(f)A<p(f) =ct. (47) 

As phase fluctuations typically occur on a much longer time- 
scale as compared to the typical gate times [44], they con- 
tribute with slow phase drifts. Therefore, we simulated N n = 
10 3 subsequent realizations of the gate where the phase drifts 
were modeled as a Wiener process with 

c =l^: [0,0.1] (48) 

where the parameter £ p determines how large the phase drifts 
are. A more detailed discussion of the process, and the moti- 
vation for the choice of the diffusion constant, can be found in 
Appendix B. 

Let us note that our analytical study (17) predicts that the 
gate is insensitive to such slow phase drifts. In this regime, we 
can set the phase to be constant A<p(f) = A<po during the gate 
interval. Accordingly, we can use the analytic expression in 
Eqs. (18) and (19) after substituting & in ^v„e iAipo , where 
A<po e [0,2ft] is a time-independent random variable. From 
these expressions, one can see that: (i) The geometric charac- 
ter of the gate, and thus its robustness to thermal noise, is not 
altered by any constant value of A<po (i-e- the fulfillment of the 
constraints (25), such that Q.\ (f g ) = 0, is independent of Acpo)- 
(ii) The effective qubit-qubit interactions (28), which deter- 
mine the gate time, jff s = L„(J^„e iA ^)( J^e iA ^)74<5„ = 
Ysn^in^n/^Sn are independent of the phase value. There- 
fore, we conclude that for strong-enough drivings and suffi- 
ciently slow phase drifts, the gate should be robust against 
phase noise. 

In Fig. 6, we present our numerical results for the resilience 
of the gate to slow drifts in the phases of the lasers. From 
these results, we can conclude that for moderate drivings, and 
phases with drifts as big as A<f> ss 0.1 ft over 10 3 gate realiza- 
tions, the gate errors still lie well-below the stringent fault- 
tolerance threshold £2 q < £ft 2 ■ We thus confirm that the pro- 
posed gate is indeed robust against a realistic phase noise, 
where the phase drifts for the time-scales of interest are well 
below A<p « 0.1ft (see the discussion in Appendix B). 




Figure 6: Phase noise: Error in the generation of the Bell state |<5~ ) 
from the initial state 1 00) ® \h\ , «i) with h n = 0, for a time evolution 
including slow noisy drifts in the laser phases. The results shown 
include a statistical average over N a = 10 3 subsequent realizations 
of the gate in the presence of the noise process (47), and we set a 
truncation of the vibrational Hilbert spaces to hold a maximum of 
"max = 7 phonons per mode. The error is shown as a function of 
the applied microwave driving strength, and for different values of 
the variance of the noise associated to the parameter £ p . Gate errors 
below the fault-tolerance threshold (FT) can be achieved £2q < £ft 2 • 

D. Partial resilience to intensity noise 

Let us now discuss the impact of possible fluctuations of 
the microwave intensity on the gate fidelities. This noise is 
modeled by substituting the microwave Rabi frequency Q.^ — > 
Q.d(t) = fid + Ai2(j(f), where An<j(f) is a stochastic process 
representing the microwave intensity fluctuations. Hence, the 
full Hamiltonian (22) should be substituted for 

n i in 

A reasonable choice for the intensity fluctuations is to con- 
sider an Ornstein-Uhlenbeck process, and set its mean and 
variance to 

ATV0=0, A£2 d (f)A£2 d (f) = f = (50) 

such that £i fixes the relative intensity fluctuations, which we 
vary in the range £j ps 10~ 4 . In addition, we assume that the 
correlation time of the intensity fluctuations is longer than the 
gate time, and set it to 1 — 1 ms. 

In Fig. 7, we integrate numerically the Hamiltonian with a 
fluctuating microwave intensity (49) and a refocusing spin- 
echo pulse at half the gate time. As seen in this figure, for in- 
tensity stabilization within the £ rj 10~ 4 regime, and interme- 
diate drivings £\i/2ft G [5, 10] MHz, the gate error is still well 
below the higher threshold Eft^ and approaches the fault- 
tolerance threshold Epj 2 for optimal drivings and the smallest 
noise £ « 0.7 • 10~ 4 . 

However, the full advantage of the strong-driving limit 
shown in Figs. 4 and 5, where £2 q < £ft 2 is attained for suf- 
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Figure 7: Intensity noise: Error in the generation of the Bell state 
|<I>~) from the initial state 1 00) ®\h\,h2) with n„ = 0, for a time 
evolution including the noise in the intensity of the microwave driv- 
ing. The results shown include a statistical average over N n = 10 3 
realizations of the noise process per point, and we set a truncation 
of the vibrational Hilbert spaces to hold a maximum of n max = 7 
phonons per mode. The error is shown as a function of the ap- 
plied microwave driving strength, and for different values of the rel- 
ative noise intensity £ £ {0.7, 1.0, 1.3} 10 4 . Gate errors below the 
fault-tolerance threshold (FT) can be achieved £2 q < EpT t for driving 
strengths Q. d /27t & 7MHz. 



ficiently strong microwave intensities, is lost due to the asso- 
ciated increase of the intensity fluctuations. Hence, it would 
be desirable to modify the scheme such that it becomes more 
robust to intensity fluctuations, while preserving its resilience 
to thermal, dephasing, and phase noise. This is accomplished 
in the following section. 



E. Doubly-driven geometric phase gates 

Inspired by recent results for prolonging single-qubit coher- 
ence times by means of continuous drivings [45], we present 
a modification of the driven single-sideband Hamiltonian (6), 
which will allow us to obtain two-qubit gates that are also ro- 
bust to intensity fluctuations of the microwave driving. We 
complement the qubit Hamiltonian (2) with a secondary mi- 
crowave driving, such that H q — > H q 2, where 



We show below that this doubly-driven model leads to a 
geometric phase gate with an additional property: it is robust 
to intensity fluctuations of the first microwave driving A£2d(f), 
while being sensitive to intensity fluctuations of the secondary 
driving AD.^(t). This already tells us that the secondary driv- 
ing should be much weaker <C Q.d, such that we can ne- 
glect the effect of its intensity fluctuations. However, in order 
to find a more detailed parameter regime for the high-fidelity 
gates, we need to explore Eq. (52) analytically. 

Based on the experience gained from the analytical study 
of the single-driving geometric phase gates (see Sec. II B), we 
move to the "dressed-state" interaction picture with respect to 
the first driving (15), which is again assumed to be on reso- 
nance with the qubit ©a = (Oq, and £2d G R. By focusing first 
on the qubit part of the Hamiltonian, which parallels the dis- 
cussion in [45], we get the following expression 

H q = L—^— a i + 

i 

+ E ( (of + iff/' cos(£V) - iaf sin(Q d ? )) e -*' + H.c. 
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(53) 



where we have introduced the detuning of the secondary driv- 
ing 5d = 0>d — (Oo- If this detuning is set correctly, the sec- 
ondary driving will decouple the qubit from the intensity fluc- 
tuations in the first line of the above equation, just as the de- 
phasing noise is minimized by the first driving (see Eq. (42) 
in Sec. IIIB). In particular, this is achieved for 8d = £ld, such 
that a rotating-wave approximation for £2d <C 4£2d leads to 

H q = i^-AQdCOof - J Ej(fid + Afid(t))of. (54) 

If we now move to the dressed-state interaction picture of the 
secondary driving, which we shall call the double dressed- 
state interaction picture 



Ud(t) = e~ if ^'' J^dOfgirE, j^dOf e irL ^0)qO- & it^„a> n ala„ 



the qubit Hamiltonian can be rewritten as 

7? q =^Ai2 d (f)(cos(in d f)cJf + sin(ii2dOffr 
-~£AQ d (;)af. 



(55) 



(56) 



H q .2 = E 5 «*>°f + 2 (iMr)cT+e- iffl i' +H d (f )<^e-^' +H.c), 

(51) 

and the secondary driving is characterized by a fluctuating 
Rabi frequency Cld(t) = Q.d+AQ.d(t ), and a frequency a>d, ful- 
filling a>d ~ (Oq, and £2d (Oq, where we have assumed again 
that £2d G K and AHd(f) is an Ornstein-Uhlenbeck process. 
This yields the doubly-driven single-sideband Hamiltonian 



#2dss = H q 2 + H v + H { 



qp. 



(52) 



where the phonon H p and qubit-phonon H qp terms correspond 
to Eqs. (1) and (4) described above. 



In this picture, it becomes clear that the noisy terms due to 
the first driving become rapidly rotating even for a weaker 
second driving provided that £i£ld <^ *C 4£2d. Hence, we 
can neglect them in a rotating-wave approximation, such that 
the qubit is only sensitive to the fluctuations of the second 
driving H q « — \ £ ; Ai2d(f)af , which have a weaker effect. 

So far, we have only treated the qubit part of the doubly- 
driven single-sideband Hamiltonian (52). The crucial point to 
address now is whether the two strong drivings can be com- 
bined with the qubit-phonon couplings responsible for the en- 
tangling gate. This question is by no means trivial, since the 
secondary driving acting in the a z -basis (54) will make the a x 
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state -dependent force (16) rotate very fast, and can thus inhibit 
the required qubit-phonon couplings. To find the correct pa- 
rameter regime, let us express the qubit-phonon Hamiltonian 
in the double dressed-state interaction picture 

i.n L 

(57) 

where we have introduced the following time-dependences 

f x (t) — cos(jCl d t) — ism.(lQ, d t)cos(Q, d t), 

f y (t) =sm(lQ d t)+icos(±Q d t)cos(Q d t), (58) 

f z (t)=sm(Q d t). 

We now study all the possibilities, (i) If the laser-induced de- 
tunings fulfill 8„ <C & d <C ^d* then all the terms in (57) be- 
come rapidly rotating, and the qubit-phonon coupling is in- 
hibited by the double driving, (ii) If 8„ « Cl d <C £l d , such 
that jl^in] <^ Cl d <^ £l d , then only the a z force can be ne- 
glected in a rotating-wave approximation. Unfortunately, the 
a x and a y forces then contribute equally, and we would lose 
the geometric character of the entangling gate (i.e. it will be 
very sensitive to thermal motion of the ions). (Hi) Finally, if 
&d < <?„ ~ ^d, such that \\&i n \ < &d, ^d, then both a*l y 
forces become rapidly rotating and can be neglected, but the 
a z force preserves its nearly-resonant character. To achieve 
this condition, instead of setting the laser frequency close to 
the bare sideband 0)l ~ (Oq — C0„ (see Sec. II A), we need to 
adjust it close to the dressed sideband resonance 

cdl^ coo-{co„-Q. d ), 8 n = a> L -{coo + £l d -ca n ). (59) 

In this regime, we obtain a single state-dependent force 

4p~-L^°f fl » e " 4 '+ H - c -> (60) 

i,n 

which will lead to a geometric phase gate in a different basis. 
We also note that the condition for the laser Rabi frequency is 
modified in this case to |£2l| -C 8(a),, — |Qd/2|). 

Let us now support the above discussion by integrating nu- 
merically the dynamics induced by the doubly-driven single- 
sideband Hamiltonian in Eqs. (57)-(58), such that the laser is 
now tuned to the regime (59). From the above discussion, we 
know that by applying a sufficiently-strong secondary driving 
£2d 3> | ^in | , we are left with a single state-dependent force in 
the o^-basis with a halved strength \^i n — > \^i n with respect 
to the previous a x -force (16). In this regime, we should re- 
cover circular phase-space trajectories when the initial state of 
the qubits is an eigenstate of the g\o\ operator. In Fig. 8, we 
analyze the strong-driving dynamics numerically. For weak 
secondary drivings [Fig. 8(a)], we observe how the different 
phase-space trajectories are not closed. Conversely, for strong 
secondary drivings [Fig. 8(b)], we recover circular orbits that 
close exactly at the expected gate time. These numerical re- 
sults support the validity of our arguments leading to Eq. (60). 



a 




Figure 8: Doubly-driven phase-space trajectories: Numerical 
analysis of the phase-space trajectories under the doubly-driven 
single-sideband Hamiltonian (57). (a) Trajectories of the forced 
normal modes for a weak secondary driving Cl d /2K = 10kHz, and 
a strong primary driving £l d /2TZ =10 MHz. Both of the phonon 
modes are in the vacuum state. In the left panel, the initial spin state 
is 1 11) or 1 00), and according to the forces (5), only the center-of- 
mass mode is driven in phase space (jc C m) = 1 /y/2(a\ +a\),(p cm ) = 
i/%/2(a] — a\). In the right panel, the initial spin state is 1 10) or|01), 
such that only the zigzag mode develops a driven trajectory in this 

case (jf zz ) = 1 /-\/2(«2 + a 2)>{Pzz) = l /^-{ a \ ~ a i)- ^ n rjotn cases > 
the trajectories are not closed, (b) Same as (a) but in the regime 
where the secondary driving is also strong Cl d /2% = 6MHz. As an 
effect of the strong driving, the trajectories are closed, and the gate 
inherits a geometric character. In the numerics we have set a trunca- 
tion of n ma x = 5 to the vibrational Hilbert spaces. 

Starting from the state-dependent force in the a z -basis (60), 
we can apply once more the Magnus expansion to get the fol- 
lowing approximation to the time-evolution operator 

U z (t e ) « ^(Oe-^W*, Hmss = l^jff^afaj, (61) 
where the coupling strengths now are given by 

J «" = -L7sc^^- (62) 

n 160„ 

Note that , in order to obtain this time evolution, it is necessary 
to neglect all other contributions coming from the Magnus ex- 
pansion. This is possible by considering the strong-driving 
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Figure 9: Entangling doubly-driven geometric phase gate: We 

display the dynamics by numerical integration of the unitary evo- 
lution induced by (57). We consider an initial state \\ffQ) = | + 
— ) ® |ni,%), where ;ii = fii = are the initial mean number of 
phonons for the cm and zz modes, and |±) = (|0) ± |l))/\/2. We 
have to set a truncation « max = 5 to the vibrational Hilbert spaces. 
In (a), we represent the dynamics of local qubit operators (of), 
which show the qubit flip | + — ) — > | — +) after t m 500/is. In (b), 
we display the fidelity between the time-evolved state and the Bell 
state I*") = (| + -) -i| - +))/v / 2- At t « 250/is, this fidelity 
approaches unity, which supports the validity of the doubly-driven 
single-sideband two-qubit gates. 



limit \ \^i n \ <C iXj <§; 4£l(ji where additional conditions simi- 
lar to the single-driving case (25) are also fulfilled 



% = r 2 ^-, 8 2 = k8 u 

01 



£l d = pd h SZ& = q8i, r,lc,p,qeZ. 

(63) 

These conditions allow us to close exactly the phase-space tra- 
jectories [Fig. 8(b)] or, equivalently, to minimize the residual 
qubit-phonon couplings that would make the gate sensitive to 
the thermal motion of the ions. We have found that these con- 
ditions can be met by imposing 



Ifld = 



215, 



Hi 



(64) 



in analogy to our older choice (32) for the single-driving gates. 
We use the same values for the detunings as in Table I with 
o"i /2 — > 5j 1 2, but now set the parameter r — 32 to optimize the 
gate fidelities, at the price of diminishing the gate speed. Note 
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that the unitary evolution (61) leads to the following table for 

t g = 7T/(87f dss ) 



(65) 



In Fig. 9(b), we represent the fidelity for the unitary genera- 
tion of the entangled state I 1 ? - ) from the initially unentangled 
state | H — ). It becomes clear from this figure, that the fideli- 
ties that can be achieved are again close to 100%. 

However, the fast and small-amplitude oscillations in the 
fidelity set an upper limit to the achieved fidelities, which is 
below the desired fault-tolerance thresholds. In order to im- 
prove the gate fidelities, we introduce a simple spin-echo refo- 
cusing pulse, such that the complete time-evolution operator 

is £/fuii(f g ,0) = t/ z ,noise(f g , | ) (o^o^ )t/ z , n oise( | , 0). In this ex- 
pression, we have also included the intensity fluctuations of 
the first driving, such that U z no i se (t2,h) is the time-evolution 
operator induced by the noisy Hamiltonian 



H 7 



£ A£l d (t ) ( cos ( \dit ) of + sin( \ Q d t ) of) 



1 

(of hit) + oJf y {t) - iof/ z (f))«„e- i5 »' +H.c, 

i,n 1 

(66) 



where all the parameters have already been introduced above. 
Note that, due to the refocusing pulse, the target entangled 
state obtained from \ -\ — ) is no longer l^ - ), but rather \^ + ). 

In Fig. 10, we represent the achieved error of this doubly- 
driven geometric phase gate for the generation of the entan- 
gled state in the presence of noise on the first driving. 
This numerical result displays the superior performance of 
this new gate, and its resilience to the intensity fluctuations of 
the first microwave driving allowing for £2 q < £ft 2 ■ At this 
point, it might be worth emphasizing the benefit of the sec- 
ondary driving. Noting that the scheme also works for clock 
states, one might ask if the impact of the intensity fluctuations 
of the microwave driving could be reduced by making faster 
gates. In fact, assuming a clock state, we can increase the 
gate speed by decreasing the parameter r from Eq. (25) with- 
out compromising the error due to the dephasing noise. In 
Fig. 11, the error for a gate with r = 2 in the presence of inten- 
sity noise Q = 0.7 • 10~ 4 on the microwave driving is shown. 
As can be seen in the figure, the error reaches E2 q ~ 1 ■ 10~ 4 , 
while the error achieved with the secondary driving is lower 
e 2q ~5-10- 5 . 

We remark that intensity fluctuations of the secondary driv- 
ing have not been considered. Such fluctuations will eventu- 
ally be the limiting factor for the gate fidelity. However, since 
the secondary driving is weaker than the first one, their ef- 
fect will be smaller. Moreover, due to the excellent fidelity 
obtained, it may also be possible to reduce the second driv- 
ing further, making the gate more insensitive to its fluctua- 
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tions. Another possibility would be to introduce a tertiary mi- 
crowave driving, which would be now detuned by the Rabi 
frequency of the secondary one (b& — ©o = This driving 
would act as a decoupling mechanism from intensity fluctua- 
tions of the second driving. We believe it should be possible 
to modify the sideband resonance (59), such that we obtain a 
triply-driven geometric phase gate in a different basis (possi- 
bly the a y basis) that enjoys a further robustness against all 
these sources of noise. Such concatenated schemes can be 
followed until the desired fidelities are achieved. 



10 Q d /2^(MHz) 30 40 

Figure 10: Intensity noise in the doubly-driven geometric phase 
gate: Error in the generation of the Bell state | > P + ) from the initial 
state | + +) ® \n\,n2) with h n = 0, for a time evolution including 
the noise in the intensity of the first microwave driving of strength 
£\l/27r = 10MHz, and relative noise of £i = 10~ 4 . We consider 
a secondary driving with £1$/2tc = 6MHz and no intensity fluctua- 
tions. The results shown include a statistical average over N n = 10 3 
realizations of the noise process, and set a truncation of the vibra- 
tional Hilbert spaces to hold a maximum of rc max = 5 phonons per 
mode. The obtained error for the doubly-driven gate is represented 
as a red star, which is well below the second threshold £2 q < £ft 2 • 
The older errors for the single microwave driving [Fig. 7] are also 
included in the background for comparison. 
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10 Q d /2^(MHz) 30 40 

Figure 1 1 : Intensity noise for a faster gate: Error in the genera- 
tion of the Bell state |4>~) from the initial state 1 00) ® |ni,«2) with 
n„ — 0, for a time evolution including the noise in the intensity of 
the microwave driving with a relative noise of § = 0.7 ■ 10~ 4 . We 
set the parameter r from eq. (25) to r = 2 which corresponds to a 
gate time t g ~ 16/xs. The results shown include a statistical average 
over N n = 10 3 realizations of the noise process, and we truncated the 
vibrational Hilbert spaces at a maximum phonon number rc max = 7 
phonons per mode. The gate error is displayed as a function of the 
applied microwave driving strength. Gate errors of £2 q ~ 1 • 10~ 4 can 
be achieved. 



IV. CONCLUSIONS AND OUTLOOK 

In this manuscript, we have demonstrated theoretically a 
two-qubit entangling gate for trapped ions, which relies only 
on a single red-sideband excitation and a strong microwave 
driving tuned to the carrier transition. By properly choosing 
the laser and microwave frequencies, we have shown analyti- 
cally that the controlled dynamics corresponds to a geometric 
phase gate in the a x basis, which we have called the driven ge- 
ometric phase gate. We have shown numerically that the gate 
is able to overcome the imperfections associated to thermal, 
dephasing, and phase noise, while achieving gate speeds com- 
parable to state-of-the-art implementations. In particular, we 
have shown that such a driven geometric phase gate can attain 
errors well below the fault-tolerance threshold gpr, ~ 10~ 4 
for sufficiently-strong microwave drivings. 

We have also analyzed numerically how this gate is still 
sensitive to intensity fluctuations of the microwave driving, 
which might be the currently limiting factor to its accuracy. To 
overcome this drawback, we have devised a new scheme with 
a weaker secondary microwave driving, which makes the gate 
robust against fluctuations in the first driving intensity. We 
have showed both analytically and numerically that, by setting 
the laser and microwave parameters in a certain regime, the 
dynamics correspond to a geometric phase gate in the a z basis. 
In this case, the limiting factor for the gate fidelity will be the 
fluctuations in the secondary microwave driving. 
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Appendix A: Magnus expansion for the driven single-sideband 
Hamiltonian 

The Magnus expansion (ME) allows us to write the time 
evolution operator of a system with a time-dependent Hamil- 
tonian H(t) as U(t,to) = where Q.(t,to) is an anti- 
hermitian operator that can be written in a perturbative series 
Q.(t,to) = Lr=i ^k{t ,to) [40]. In this work, we are interested 
in the ME to second order, where the first- and second-order 
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terms are given by 



fli(f,f ) = -i / dfifffa) 

J to 

£l 2 {t,t ) = [ dti f 1 6t 2 [H(h),H{t 2 )]. 



(Al) 
(A2) 

In Sec. II B, we used the results for the second-order ME for 
the qubit-phonon Hamiltonian (14), namely 

jr. 

// qp =^^(aj : +icr]'cos(£2 d r)-iff;"sin(£2 d f))fl„e" i6 " r +H.c., 

(A3) 

which is expressed in an interaction picture with respect to the 
microwave driving. In this Appendix, we present a detailed 
derivation of the ME and discuss its more relevant terms. 

Let us start by considering the first-order contribution, 
which is calculated directly from integrals of the qubit-phonon 
Hamiltonian, and leads to the following result 



fli(f,0)=£ 

,i(O d -S„)r 



2 



-iS„t _ 



I ) -^-<JjCL n 



H.c., 
(A4) 



which coincides with Eq. (18) in the main text, and contains 
the different state-dependent forces. The second-order term 
can be split into three different parts, namely 

Q.2 (r , 0) = Q. a 2 (t , 0) + £l\ (r , 0) + Q§ (r , 0) . ( A5) 

The first part Q%(t , 0) contains the terms that are linear in time, 
and can be interpreted as a spin-spin Hamiltonian generated 
by the three non-commuting state-dependent forces 



(A6) 

where 5» is the Kronecker delta and the coupling constants 
Tjf, M*jf, ICjf and the driving corrections AQ.j„ are given by 
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n d — s„ 


i2 d 4 
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f n d 


+ 8 n/ 



(A7) 



At this point, let us note that our specific configuration con- 
siders a propagation direction of the laser beams fulfilling 



kL ■ r j = 0, such that the forces J£ )„ are purely imaginary 
(see Eq. (5)), and the couplings M*jf vanish exactly. More- 
over, in the strong-driving limit 8 n <C £2 d , we find that K^f °= 
^jn^lJD^, which can be neglected to order ^(|), with 
4 = (ilL^Jn) 2 /^ *C 1, as was done in Eq. (19) of the main 
text. Altogether, we conclude that the most important part is 



E n A%„ffj%)' (A8) 



Let us now move to the second part il| (t , 0) , which contains 
oscillatory couplings between the internal qubit states 



yeff 

Q|(t,0) = iE 4^- (sin(iV)crj + (1 - cos(£2 d f)) aj 



f4n d 



fe- i5 "'-l) 



" ? I ~ T" V a >*" 

j,k,n L 

+^ n (e'^-^)'-l)+4„(e-'^+^)'-l) 



H.c. 



(A9) 



where we have introduced the two-body spin operators 



jkn 



jkn 



1 * v 



1 



2<5„ (£2^ - 
1 



25„(n d -<5„) 



{-iojol-ojo*), 



1 
1 

2<5„(i2 d 



(A10) 



Let us now discuss which are the leading terms in the above 
expression (A9). The first line of Eq. (A9) only contributes 
perturbatively since J*jf -C 8fi d . From the remaining contri- 
butions to (A9), note that most of the terms in (A10) scale 
with ■J^j„JP r £ n /Q.d8„, and can be thus neglected to order 
where % = (£2lJ7„) 2 / £2 d 5„ -C 1, as was done in Eq. (19) of the 
main text. Finally, considering again that the forces are purely 
real in our setup, we get the following leading contribution 



nfoo)«-i£ 

j,k 



^ jn ^* n sm(8 n t) 

452 



<jj<r k . 



(All) 



The third and final part Q|W contains contains the residual 
spin-phonon couplings, and reads as follows 

Cl c 2 (t)=itl t ACl jn a]ata n + £ {f1f{t)ala n -H.c.) 

j,n j,m^n 



J.m,n 



(A12) 



where A£lj n was introduced in (A7), and we have introduced 
the following time-dependent single-qubit operators 



1 



1 



8(5„ - 5 m ) \ £2d - 8m + 8„ 



-i(<5„-5,„)f _ 



■Dai, 



(A13) 
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together with 

fir® = b) mn (e- [s »> - 1) +b) mn ^- i{SM 1) 

+b) mn {e^- s »- 5 ^ - 1) +b% m (e-^ +s » +5 ^ - 1), 

(A14) 

and, finally, 

hf (0 = c jmn ( 



e -i5„r _ j-j 



-4„„(e^- 5 »)'-l)+4„„(e-'^ +5 »)'-l) 



+ 4„„(e i ( n <' +5 <»- s ")' _ l) +4„„(e- i ( a <' +5 "- 6 '"> - 1). 

(A15) 

In these last two expressions, for notational convenience, we 
have introduced the following list of single-qubit operators 
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(A16) 



s?- 3P* 

,sr jn^ jm 



85„(£2 d 



(A17) 

Once all the expressions of the Magnus expansion have been 
described, let us analyze the leading-order contribution to 
the dynamics. First of all, note that in the strong-driving 
limit i2 d 3> 8 n , all the expressions in Eqs. (A16) and (A17) 
scale as & j n &*- m ] £l$8 n , which is again on the order of 

where % = {Q.\J] n ) /i2 d 5„ <C 1, and can be directly neglected. 
Hence, the leading contribution to this part is 

Q.\{t) « k£AO yn aj=4a n + £ C/5?(0«U - Hx.). 

(A18) 



It is now easy to convince oneself that all the leading contri- 
butions to the second-order (A8), (All) and (A18), form the 
expression (19) used in the main text. 



Appendix B: Stochastic processes for the noise sources 

In this part of the appendix, we present a more detailed dis- 
cussion of the origin of the dephasing, phase, and intensity 
noise in ion trap setups. As described in the main text, all 
these noise sources can be modeled by a particular stochastic 
Hamiltonian. We use this appendix to introduce the particular 
stochastic models used in this work, and their main properties. 

We use two continuous memoryless stochastic processes, 
also known as Markov processes, namely the Ornstein- 
Uhlenbeck and the Wiener process. In the following, we intro- 
duce briefly both processes, and present some of their proper- 
ties that will be important for our numerical simulations. Sub- 
sequently, we will describe how they are used to model the 
fluctuations of a specific quantity in the trapped-ion setup. 

a) Ornstein-Uhlenbeck process.- The Ornstein-Uhlenbeck 
(O-U) process (9(f) is determined by a diffusion constant c and 
a correlation time T, and evolves according to the following 
stochastic differential equation 



dO(t) 
dt 



1 



0{t) + VcT(t) 



(Bl) 



where T(t) is a Gaussian white noise. Remarkably, equa- 
tion (Bl) is exactly solvable [41], which allows to show that 
the O-U process is a normal process with mean and variance 



-t/t 



Var{0(f)}=^(l-e- 2 '/ T ). (B2) 



0(0 = O(0)e 



As one can see in (B2), the correlation time T sets the time 
scale on which the asymptotic values of the process are 
reached. Assuming a zero mean, the auto-covariance of the 
O-U process is given by 



cz 



O(0O(0) = ye-'/*(l- 



-2//T 



(B3) 



where one can see that z also governs the time-scale over 
which the noise is correlated. Finally, the essential ingredi- 
ent for the numerical simulations is the exact update formula 



0(t+At) =0(t)e 



-Af/T_ 



cz 



(1-e 



-2Ar/T\ 



1/2 



(B4) 



where n is a unit Gaussian random variable. We note that this 
formula is valid for any finite time step At, which will turn out 
to be very useful for the numerical simulations. 

b) Wiener process.- In the limit of very large correlation 
times z — > oo, an O-U process with diffusion constant c be- 
comes another type of Markov process, which is known in the 
literature as a driftless Wiener process W(t) [41]. The stan- 
dard form of the Langevin equation for the Wiener process 
W(t) reads as follows 



W(t + Af) = W(t) + (cAt) l/2 r 



(B5) 
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where n is again a unit Gaussian random variable. Note that 
this serves again as an update formula that will be very use- 
ful for our numerical simulations. Introducing once more the 
Gaussian white noise function T(t), equation (B5) can be re- 
cast as a stochastic differential equation 



dW(t) 
dt 



VcT(t), 



(B6) 



which corresponds formally to Eq. (B 1 ) in the limit T — > °°. 
The Wiener process is also a Gaussian process and is charac- 
terized by a mean and variance 



W(t)=W(0), Var{W{t)} = ct. 



(B7) 



In the following paragraphs, we will discuss the concrete 
application of the above stochastic processes for the simula- 
tions of dephasing, phase, and intensity noise. For the numer- 
ical simulations, the values of 0(t), W(t) are sampled accord- 
ing to the update formulas in Eqs. (B4) and (B5). Note that 
for all simulations we set (9(0) = W(Q) = for the stochastic 
processes. 

i) Dephasing noise - The internal structure of the atomic 
ions can be perturbed by uncontrolled electric or magnetic 
fields. Magnetic-field fluctuations are typically leading, and 
are expected to be the major source of dephasing [11], which 
ultimately leads to the loss of coherence. To understand the 
effect of a fluctuating magnetic field, note that the Zeeman 
shift changes the qubit resonance frequency by an amount 
Ao)o given by 

A(0 = d B (o \ Bo (B-Bo) + ^dl(Oo\ Bo {B-B ) 2 , (B8) 

where the field B fluctuates around the average value Bo- 
Many hyperfine transitions are first-order field insensitive at 
B = Bq = 0, but the need to resolve the hyperfine structure 
for the definition of our qubit (see Fig. 1) requires a finite 
external magnetic field. Interestingly, there is a family of 
states/transitions, often referred to as clock states [46], where 
the linear Zeeman shift vanishes at a certain value Bo, yielding 
thus a higher resilience to magnetic -field noise. 

Although our scheme applies equally well to such clock 
states, we would like to consider magnetic -field sensitive 
states, and show that the strong driving also protects the qubits 
from magnetic-field noise. Therefore, we will assume that 
Aa)o(f) = ^ ^MfiB(f), where g is the hyperfine g-factor, jUg 
the Bohr magneton, and B(t) a randomly fluctuating magnetic 
field. Thus, Aoo(f) also fluctuates randomly, and the noise is 
described by a stochastic Hamiltonian term 



//flu 



(B9) 



as introduced in Eq. (38) of the main text. For clock states, 
Aa)o(f) will contain the weaker quadratic Zeeman shift, and 
the ac-Stark shifts due to fluctuating laser intensities. 

In order to reproduce the exponential decay of the coher- 
ences typically observed in experiments, we modeled the fluc- 
tuations of the resonance frequency Acoo(f) as an Ornstein- 
Uhlenbeck (O-U) process (Bl). According to Eqs. (B2) the 
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Figure 12: Dephasing noise process: Fluctuations of the qubit reso- 
nance frequency Acoo(f) during a single realization of the driven geo- 
metric phase gate t g ~ 63 /is. The dashed lines represent the standard 
deviation that follows from (B10) of the O-U process that models the 
dephasing, while the solid line represents the stochastic process it- 
self. We set the decoherence time to T 2 = 25 /is, which fixes the 
correlation time T = 2.5 /is, and the diffusion constant c = 2/(T2T 2 ) 
of the O-U process. 



fluctuations in the resonance frequency are characterized by 



CT 



Aa)o(f)=0, Var{Aa>o(f)} = Aa>o(f) : 

(B10) 

(recall that we set Ao(0) = 0) and following (B4) the update 
formula reads 



A(0o(f +Af) = Ao>o(0 e 



-At ft _ 



CT 



-a 



-2A?/t 



1/2 



(Bll) 

In our noise model, the temporal decay of the qubit coher- 
ences is given by 

(o?(t f )) = (a*(0))e- 1 2( Y2 ( t f)), (B12) 

where T(ff) is the integral of the O-U process [41]. The exact 
expression for (Y 2 (tf)) is given by 



•t 4 

(Y 2 (t f )) = J dt'/ dt"(Acoo(t')Ac%(t")) 





-CT 2 



(?f-2T(l 



-f f /t\ 



|(l-e- 2 ^)) 



(B13) 



In the limit tf ^> T, and using Eq. (B13), we get (af(tf)) = 
(of(0))e~' f,/r2 , where we have introduced the dephasing time 



T 2 



CT 



2 ■ 



(B14) 



For all numerical simulations, we set T = O.I/2. Thus, af- 
ter choosing T2 the positive constant c is determined by Eq. 
(B14). As we shall show in the following, it is required that 
the time step At used in the numerical simulations be much 
smaller than T. 
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Figure 13: Phase noise process: Fluctuations of the laser phases 
A<p(/) leading to a slow drift during subsequent repetitions of the 
driven geometric phase gate, ff = 10 3 / e ~ 63 ms. The dashed lines 
represent the standard deviation that follows from (B21) of the 
Wiener process that models the phase drifts, while the solid line rep- 
resents the stochastic process itself. We set the relative phase fluc- 
tuations to be £ p = 1CP 2 , which fixes the diffusion constant of the 
process according to Eq. (B22). 
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Figure 14: Intensity noise process: Fluctuations of the microwave 
driving intensity A£l^(t) during a single realization of the driven 
geometric phase gate t„ ~ 63 jls. The dashed lines represent the 
standard deviation that follows from (B2) of the O-U process that 
models the intensity fluctuations, while the solid line represents the 
stochastic process itself. We set the relative intensity fluctuations to 
be £i = 1.3 • 10~ 4 for a driving of Q.^/2n ~ 7MHz, the correlation 
time was fixed to z = 1 ms, which fixes the diffusion constant of the 
process according to Eq. (B26). 



Since time is discretized in the numerical simulations, the 
integral is approximated by a sum 



(Y 2 (k)) 



CX 



Ate t > Are * 

n=l n'=\ 
CT N nU N 

^Afe 2 ^ Y_, Afe~ 



(B15) 



n=l n'=n+l 

where NAt = ff. After some algebra, equation (B 15) gives 



(Y 2 (tf)) 



CT 
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(B16) 



1 -cosh(Af/r) 



In the limits At -C T, and NAt = tf^$>X, we can approximate 
Eq. (B16) by 



(y 2 (f f )) w ct 2 (NAt - t) « ct\ 



(B17) 



giving a good approximation to the analytical result in 
Eq. (B14). A similar result is obtained for At fa X. On the 
other hand, for the case Af 3> x, we get 



(F 2 (f f )> a °?-Att f 



(B18) 



such that the exponential decay of the coherences will be 
much faster than expected. Consequently, we chose the time 
step Af in the numerical simulations to be Af -C x. 

In Fig. 12, we represent the fluctuations of the qubit reso- 
nance frequency given by an O-U process fulfilling the above 
requirements. The time-scale in this figure corresponds to a 
single realization of the driven geometric phase gate, and we 
typically average over A^ n = 10 3 gate realizations. 



ii) Phase noise.- Another source of noise in trapped-ion 
experiments is due to the fact that the phases of the lasers 
at the positions of the ions are not constant, but rather sub- 
jected to slow drifts [47]. To study the effects of phase fluctu- 
ations on the driven geometric phase gate, we make the sub- 
stitution J£",„ — > ^ ^ ■ n e lA(^ '' ^ • , for the laser-induced sideband cou- 
plings such that the qubit-phonon couplings (3) become 



(B19) 



as introduced in Eq. (46) of the main text. As reported in ex- 
periments [44], a phase shift of 2n takes place on time-scales 
on the order of ~ 10s. Since these drifts occur on a very long 
time-scale, it appears convenient to model these fluctuations 
as a Wiener process. Therefore, the phase fluctuations evolve 
according to Eq. (B5) as 



Acp(t+At) = A<p(f) + (cAf) 1/2 n. 
The mean and the variance of the process are given by 



(B20) 



A<p(r) =A<p(0) =0, Var{A<p(f)} = A<p 2 (f) =ct, (B21) 

In the numerical simulations, we did not rely on the as- 
sumption that the phase of the lasers is constant during one 
gate realization (see the qualitative discussion in the main 
text). Since the fluctuations during one gate realization are 
very small, the relevant question is if the phase fluctuations 
affect the gate performance for a large number of consecutive 
gates. Therefore, we simulated A^ n = 10 3 consecutive gate re- 
alizations with fluctuating sideband couplings ^j„e l ^'\ 

The stochastic process A<p(f) describing the phase fluctua- 
tions is then fully characterized by its diffusion constant c and 
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its initial value A<p(0) = 0. With the results of [44], and an ex- 
pected gate time f g m 63 fls, we can estimate c following Eqs. 
(B7)as 



(C P ^) 2 

10 3 r e 



C P «o.oi. 



(B22) 



The noise process A<p(f) can then be generated using the up- 
date formula (B20), which demands the time step At for the 
simulation to be sufficiently small. This requirement should 
be fulfilled for our particular choice ofN = 200 steps per gate. 
In Fig. 13, we represent the phase fluctuations modeled by a 
Wiener process fulfilling the above requirements. 

Hi) Intensity fluctuations - The scheme for driven geomet- 
ric phase gates introduced in this paper achieves very high 
fidelities for entangled two-qubit states in the regime of very 
high powers of the microwave driving. Unfortunately, in this 
regime, the output amplitude of microwave sources cannot be 
held perfectly constant, and we should consider fluctuations 
in the microwave intensity. This leads to a noisy qubit Hamil- 
tonian 

N 

// q = ^3 c °o (J "+2(( £2 d + An d (0)c7+e- iffl d'+H.c.), (B23) 



(recall that we set Ai2 d (0) 
mula 



0) and the following update for- 



AH d (f +Af) = An d (f)e 



-Af/T 



CT 



(1 



-2M/-C 



1/2 



(B25) 

We expect that typical correlation times for such microwave 
intensity fluctuations will be much larger than the expected 
gate times. Accordingly, z = 1 ms is set to be much larger than 
the gate time f g = 63 /is. Moreover, we have assumed relative 
intensity fluctuations on the order of £j = A£2 d /i2 d « 10~ 4 , 
which requires a very accurate stabilization of the microwave 
sources, but is in principle possible (cite). Identifying the rela- 
tive fluctuations of the microwave with the standard deviation 
of the O-U process, the diffusion constant c of the O-U pro- 
cess is given by 



2(CAi) 2 r 
c= , Qi' 



10" 



(B26) 



where AH d (f) represents the fluctuations of the microwave 
Rabi frequencies introduced in Eq. (51) of the main text. We 
model these fluctuations by an O-U process (B 1 ), character- 
ized by 



AQ d (r)=0, Var{AH d (f)}: 



2~ 



-2t/T 



(B24) 



In Fig. 14, we represent the fluctuations of the microwave 
Rabi frequency given by an O-U process fulfilling the above 
requirements. The time-scale in this figure corresponds to a 
single realization of the driven geometric phase gate, and we 
typically average over N n — 10 3 gate realizations. 
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